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ABSTRACT 

We present high-quality Very Large Array (VLA) images of the Fanaroff & Riley Class 
I radio galaxy 3C31 in the frequency range 1365 to 8440 MHz with angular resolutions 
^ ' from 0.25 to 40 arcsec. Our new images reveal complex, well resolved filamentary 

substructure in the radio jets and tails. We also use these images to explore the spectral 
structure of 3C 31 on large and small scales. 

We infer the apparent magnetic field structure by correcting for Faraday rotation. 
Some of the intensity substructure in the jets is clearly related to structure in their ap- 
parent magnetic field: there are arcs of emission where the degree of linear polarization 
increases, with the apparent magnetic field parallel to the ridges of the arcs. 
. The spectra of the jets between 1365 and 8440 MHz are consistent with power 

laws within 60 arcsec of the nucleus. The spectral indices, a (flux density oc v~ a ) are 
significantly steeper (a = 0.62) within «7 arcsec of the nucleus than between 7 and 
50 arcsec (a = 0.52 - 0.57). The spectra of the arcs and of the jet edges are also slightly 
flatter than the average for their surroundings. At larger distances, the jets are clearly 
delimited from surrounding larger-scale emission both by their flatter radio spectra 
and by sharp brightness gradients. 

The spectral index of 0.62 in the first 7 arcsec of 3C31's jets is very close to that 
, found in other FRI galaxies where their jets first brighten in the radio and where 

' X-ray synchrotron emission is most prominent. Farther from the nucleus, where the 

spectra flatten, X-ray emission is fainter relative to the radio. The brightest X-ray 
emission from FRI jets is therefore not associated with the flattest radio spectra, but 
with a particle-acceleration process whose characteristic energy index is 2a + 1 = 2.24. 
The spectral flattening with distance from the nucleus occurs where our relativistic jet 
models require deceleration, and the flatter-spectra at the jet edges may be associated 
with transverse velocity shear. 

Key words: galaxies: jets - radio continuum: galaxies - magnetic fields - polarization 
- MHD - acceleration of particles 



o 



> 
On 

in 



m 
O 
oo 
o 



1 INTRODUCTION 



The radio galax y 3C 31 is a nearby m ember of the low- 
luminosity FRI l|Fanaroff fc Rilevl 1 1974 ) class. Its promi- 
nent, symmetrical twin radio jets were disc overed bv | Burch 



1977) and imaged at higher resolution by iFomalont et al 
1980). The jets are asymmetric on scales up to «10 kpc and 
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only the brighter jet is detected on parse c scales |Lara et al.l 
I l997fl More diffuse tails of emission (|Strom et all 1 19831 : 
iJaeers! Il987l ) extend for ~20 arcmin on either side of 
the nucleus. The source shows strong linear polarization 
and a pronounced depolarization asymmetry (lBurch|["l979l ; 
IStrom et alii 19831 ). 

The d etection of the brighter jet a t opti cal wave- 
lengths by iButcher. van Breugel fc Milevl (Il980t) was not 
confirmed by Fraix-Burnet. Golombek fc Macchettol (|l99ll ) 
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or lKeell (|l988l ), but lCroston et all (|2003l ) presented evidence 
for associated optical emission between 5 and 10 arcsec from 
the nucleus. X-ray em ission from this jet was found by 
IHardcastle et all (|2002h . 

3C 31 is identified with the bright elliptical 
galaxy NGC 383, which we take t o have a r e dshift 
of 0.0169 (the mean value from Smith et al. _ 200(1 



Huchra. Vogelev fc Gellerl Il999l and [5e Vaucouleurs et al.l 
19911 ). N GC 383 is the brightest member o f a rich group of 



galaxies (jArdl 19661 ; IZwickv fc Kowallll96Sl ) and, in common 
with other nearby radio galaxies, has a dusty nucleus 
ijMartel et al.lll999h . Hot gas on gro up and galactic scales 
has been detected by X-ray imaging (|Komossa fc Bohringerl 
1 19991 ; IHardcastle et alJliooa ). 

Deep 8440-MHz Very Large Array (VLA) images of 
3C 31 at resolu t ions of 0.75 and 0.25 arcsec were made by 
ILaing fc Bridle! l|2002al ). They showed that the total inten- 
sity and linear polarization within 30 arcsec of the nucleus 
can be modelled on the assumption that the jets are in- 
trinsically symmetrical, axisymmetric, relativistic and de- 
celerating. They also assumed that the jet flows are station- 
ary, in the sense that they are characterized by the same 
global parameters at all times. By fitting parameterized 
models of the velocity field, emissivity and field-ordering, 
they were able to deduce that the angle to the line of sight 
is ~52°, to quantify the deceleration and to establish that 
transverse velo c ity gra dients must be present in the jets. 
ILaing fc Bridie] l|2002bl) used this kinematic model and a 
description of the exte rnal gas pressure and density from 
IHardcastle et all (2002) to make a conservation-law analy- 
sis. They showed that the inferred velocity field is consistent 
with jet deceleration by entrainment of thermal matter and 
deduced the energy flux and the variations of pressure, den- 
sity and entrainment rate with distance from the nucleus. 
ILaing fc Bridle! l|2004l ) examined adiabatic models for the 
jets in 3C 31: these models predict too steep a brightness de- 
cline along the jets for plausible variations of the jet velocity 
and fail to reproduce the observed magnetic-field structures. 

In this paper we describe VLA observations of 3C 31 ob- 
tained under NRAO observing proposals AF236 and AL405 
at frequencies between 1365 and 8440 MHz and resolutions 
of 0.25 - 40 arcsec (the 8440-MHz observations are those 
presented by ILaing fc Bridle]|2002al ). We then discuss three 
aspects of the astrophysics of 3C31 whi ch complement 
the tre at ment of its inner jet dynamics by ILaing fc Bridla 
|2002al lbl. l2004h as follows. 



well-defined outer e dges). As in other members of the class 
|Parma et al.l|l9 99), its r adio spectrum steepens with dis- 
tance from the nucleu s l|Burchl 1 19771 ; IStrom et all Il983l ; 
lAndernach et all Il992l ). Analyses of such spectral varia- 
tio ns from low-resolution o bservations (as made for 3C 31 
by lAndernach et all Il992h interpreted them as the re- 
sults of radiative losses on an energy spectrum initially 
of power-law form and derived corresponding time-scales, 
which were identified with d ynamical ages. More recently , 
high-resolution observ ations (Kat z-Stone fc Rudnickl 1 19971 ; 
IKatz-Stone et al.ll 19991 ) have shown that spectral structures 
of a number of FRI sources are better described as a flat- 
spectrum 'spine' surrounded by a steeper-spectrum 'sheath' 
rather than as a single component whose spectrum varies 
smoothly with distance from the nucleus. In 3C31, we have 
good spatial resolution over a range of frequencies and can 
examine the spectral variations in more detail. 

(iii) Jet base spectra, X-ray emission and particle acceler- 
ation processes. We have recently shown that accurate imag- 
ing of radio spectra in FRI jet bases can reveal correlations 
between spectral index, X-ray emission and t he jet kine- 
matics as inferred from our relativistic models l|Laing et al.l 
l2006al lbl). The shortness of the synchrotron lifetimes for elec- 
trons radiating in the X-ray band requires them to be accel- 
erated in situ, so there must be a direct connection between 
the radio spectral index and the particle acceleration pro- 
cess. For 3C 31 we can determine the radio spectrum very 
accurately at high resolution and examine the connection 
with X-ray emission and kinematics in greater detail. 

The observations and their reduction are described in 
Section[2]and the images are presented in Section[3] Discrete 
features in total intensity and the apparent magnetic-field 
structure are discussed in Section [4] and variations in radio 
spectrum across the source are described in Section [S] Sec- 
tion [6] summarizes our main conclusions. Faraday rotation 
and depolarization of t he polarized emiss ion are analysed 
m a companion paper (|Laing et al.1 12008) which discusses 
the structure of the magnetic field in the foreground screen 
that is responsible for the Faraday rotation measure (RM) 
fluctuations. 

We assume a concordance cosmology with a Hubble 
Constant H = 70 km s" 1 Mpc -1 , Q A = 0.7 and tt M = 0.3. 
At the redshift of NGC 383 (z = 0.0169), this gives a linear 
scale of 0.344 kpc/arcsec. We also adopt the sign convention 
S(v) oc v~ a for the spectral index, a. 



(i) Arcs and filaments. The jets in 3C 31 contain discrete 
features - arcs and non-axisymmetric knots - visible in both 
total and linearly polarized intensity. Such features, and fil- 
amentation in general, are common in extragalactic syn- 
chrotron sources when observed with sufficiently high res- 
olution and good sp atial frequency coverage. In particular, 
ILaing et~al . (2006b) found similar arc-like structure in the 
twin-jet source 3C 296 and suggested that systematic differ- 
ences between the arcs in the main (brighter) and counter- 
jets could also be a manifestation of differential relativistic 
aberration. In 3C 31, we study how the intensity gradients in 
such features are related to their magnetic-field structures. 

(ii) Large-scale spectral structure. 3C 31 belongs to the 
sub-class of FRI sources with tails of diffuse emission ex- 
tending away from the nucleus (as opposed to lobes with 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Observations 

The observations were made with the NRAO VLA at fre- 
quencies between 1365 and 8440 MHz and in all four config- 
urations of the array. A journal of observations is given in 
Table [1] The flux-density scale was set using observations 
of 3C 48 or 3C 286 and the zero-point of E-vector position 
angle with reference to 3C 138 or 3C 286. A single pointing 
centre (coincident with the galaxy nucleus) was used for all 
observations. 
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Figure 1. Low-resolution, mean L-band images (Section I2.2[l optimized to show the extent of the observed large-scale structure, (a) 
Contours at 40arcsec FWHM. The levels are -1, 1, 2, 4, 8, 16, 32, 64, 128 X 3.5 mjy beam" 1 , (b) Grey-scale at 5.5arcsec FWHM. The 
displayed range is — lmjy beam -1 . Key regions are labelled, as discussed in the text. Note that the primary beam correction leads to 
increased noise at the edges of the field and that the low-level stripes in the tails (panel b) are artefacts. Chromatic aberration is large 
at the edges of panel (b), so images of point-like sources are smeared radially by a factor rj3 and reduced in peak flux density by the 
same amount. The effect on the brightness distribution of 3C31, which is fully resolved on these scales, is slight (Section 12.21 1. 



2.2 Data reduction 

Calibration and imaging of the VLA data were done us- 
ing the AlPS software package and followed standard proce- 
dures with o ne exception: we used the iterative technique 
described bv lLaing fc Bridle! l|2002ah in order to correct for 
core variations and amplitude errors when combining data 
from different configurations at 4985 and 8440 MHz. For the 
lower-frequency observations, core variability was below the 
level we could measure and no corrections were made. 

We used both clean and maximum entropy deconvo- 
lution methods for the total-intensity images and compared 
the results (in the latter case, we subtracted the core before 
deconvolution and restored it again afterwards). In all cases, 
the differences between the images (after convolution with 
the same Gaussian restoring beam) were less than the larger 
of 3cr_f (Table [2} and 0.01/. The mean levels were consistent 



to <C oi and the total flux densities differed by <1%. As 
usual, the maximum entropy images of extended, uniform 
brightness distributions were locally smoother and showed 
little of the patchy structure or high-frequency ripples char- 
acteristic of the CLEAN algorithm. We therefore show the 
maximum-entropy 7 images at resolutions of <1.5arcsec. 
All quantitative measurements have also been checked on 
the CLEAN images. 

The parameters of the images used in this paper are 
summarized in Table [2] All of the images have been cor- 
rected for the effects of primary beam attenuation: as a con- 
sequence, the quoted rms noise levels strictly refer to off- 
source regions close to the pointing centre. We have used 
the mean of the images at 1365, 1435, 1485 and 1636 MHz 
to provide several of the illustrations, and we refer to these 
as 'mean L-band images' (using the radio vernacular for des- 
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Figure 2. Grey-scale representations of the mean L-band image of 3C31 (Section I2.2jl at a resolution of 5.5 arcsec FWHM. (a) Total 
intensity. The grey-scale range runs from — 8 mjy beam -1 , (b) A Sobcl-filtered version of the same image, normalized by the total 
intensity (see text). 



ignation of this frequency band). We have checked the regis- 
tration of the images at different frequencies using the peak 
of the core emission and isolated point sources in the field. 
We believe that the images at f .5 and 5.5 arcsec FWHM are 
aligned to better than FWHM/20. 

The effects of limited short-spacing coverage and pri- 
mary beam attenuation are important at the two high- 
est frequencies, particularly for spectral analysis. We esti- 
mate that the 4985-MHz and 8440-MHz images are ade- 
quately sampled for scales < 300 arcsec and < 140 arcsec, 
respectively, i.e. ±150 and ±70 arcsec from the nucleus 
l|Tavlor. Ulvestad fc Perievl 120031 ). There is less power on 
large spatial scales in Q and U, and we estimate that these 
are reliable at significantly larger distances from the nucleus. 

Chromatic aberration ("bandwidth smearing") causes 
point-like sources far from the phase centre to appear dis- 
torted and reduced in peak flux density. This effect is large at 
the periphery of our 5.5-arcsec L-band images (Fig.[]J>) and 
typically causes a loss of peak intensity of 5 - 10% for point 
sources at the edges of the areas we sho w in the remaining 
figures l|Tavlor. Ulvestad fc Perley 2003). Chromatic aber- 



ration preserves flux density, however, and 3C 31 is fully re- 
solved on all scales where it is significant. Its effect on the 
images of the source shown in this paper is therefore very 
small compared with other errors. A 10-s integration time 
was used for all of the observation s, so time-average smearing 
is negligible over all of our fields (|Tavlor. Ulvestad fc Perievl 
l2003h . 

When deriving images of polarized power, P — (Q 2 + 
U 2 ) 1 ^ 2 , and degree of polarizatio n p = P/I, we made a first - 
order correction for Ricean bias (|Wardle fc Kronberd[l973 ). 



3 THE IMAGES 

3.1 Large-scale structure 

Fig.[T]shows the large-scale structure of 3C 31. The 40-arcsec 
FWHM image (Fig. [1J,) shows the maximum detectable ex- 
tent of the structure at high signal/noise ratio. The low- 
est contour is essenti ally identical to that seen in previous 
low-frequency images (lBurchlll977l ; IStrom et al 1ll983l ; |jagersl 



1 19871 ; lAndernach et al] Tl992) 
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Figure 3. Grey-scale images of 3C31 at a resolution of 1.5 arcsec FWHM. (a) Mean L-band maximum-entropy deconvolution (Section l2.2l l, 
— 1 mjy beam -1 , (b) 4985 MHz maximum entropy deconvolution, — 1.5 mjy beam -1 . The area covered by (b) is shown by the box 
on (a). 



In Fig. Hlb), we show a grey-scale representation of the 
mean L-band image at a resolution of 5.5 arcsec FWHM, 
over a flux-density range chosen to emphasise the low- 
brightness structure. All of the emission visible at lower res- 
olution is seen on this image, albeit at a lower level of signifi- 
cance and with some artefacts. Key regions of the source are 
labelled: we refer to North and South jets, spurs and tails. 
The essential distinction between tails and spurs is that the 
former extend away from the nucleus whereas the latter ex- 
tend back towards it, at least in projection. Filamentation 
and fine structure are already apparent in Fig. [T^b) : the bi- 
furcation of the North tail and the filaments in the North 
spur are best seen on this display. The remarkably straight 
East edge of the South spur is more prominent than on ear- 
lier images. 

In a display of the mean L-band image at 5.5 arcsec 
FWHM resolution over a larger intensity range (Fig. [2^,), 
the true extent of the jets becomes apparent. In order to 
display the intensity gradients more clearly, we show a Sobel- 
filtered r epresentatio n of the image in Fig. [2jb) . The Sobel 
operator (|Prattlll99ll ) computes VI | and therefore enhances 
brightness gradients. We show an image of fractional inten- 
sity gradient |V7|/J. This image shows that the boundaries 



of the jets are marked by strong fractional intensity gradi- 
ents. These can be followed to the positions marked N and S 
in Fig. [2jb) and we therefore use them to define the bound- 
aries of the jet, spur and tail regions. At the points N and S, 
there is evidence for curved features with sharp brightness 
gradients slightly inside the outer boundary of the source. It 
is tempting to associate these features with the redirection 
of the jets to form the tails and spurs. It may also be that 
the spur and tail are parts of the same physical structure, 
seen in projection on either side of the jet. 



The South jet has a well-collimated region of compar- 
atively low emissivity (labelled 'A' in Fig. |2ji) which be- 
gins where it becomes superposed on the spur. Thereafter 
it bends through m 90° and brightens before widening and 
appearing to merge with the spur and tail emission. The 
North jet makes a sharp double bend before widening. Both 
jets show complex, but well-resolved sub-structure after they 
widen. This sub-structure is clearly evident in Fig.(2ja), but 
is not prominent in the fractional intensity gradient image 
(Fig.(2}j). The straight edge in the South spur noted earlier is 
particularly prominent on the Sobel- filtered image (Fig. [2)3). 
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Figure 4. Grey-scale image of 3C31 observed at 8440 MHz with a resolution of 0.75 arcsec FWHM. (a) Total intensity (0-1 mjy 
beam -1 ), (b) Sobel-filtered and normalized by total intensity. The 'arcs' discussed in Section 13.31 are indicated and labelled with their 
morphological types where these can be ascertained. 



3.2 Total intensity images at 1.5 arcsec resolution 

A clearer representation of the termination of the jets is 
given by the mean L-band image at 1.5-arcsec FWHM reso- 
lution (Fig.|3^.), in which the tail emission is heavily resolved. 
The internal structure of the jets starts to be revealed by the 
4985-MHz image at the same resolution (Fig. [3j>) . The main 
new features are as follows: 



(i) There is evidence for considerable filamentation and 
other sub-structure within both jets. We discuss this in de- 
tail below (Sections 13.31 and 14. 2[) . 

(ii) The low-emissivity region A of the South jet men- 
tioned earlier is seen in more detail. It appears to narrow 
slightly before bending and brightening. 

(iii) The region where the South jet expands inside the 
South tail shows complex, resolved filamentary fine struc- 
ture. 



3.3 8440 MHz total intensity images at 0.75 and 
0.25 arcsec resolution 

The higher-resolution i mages shown here cover a wider area 
than those discussed by lLaing fc Bridle! (1200281 ). The 0.75- 
arcsec FWHM image (Fig. [4^,) shows a number of features 
with sharp intensity gradients. Examples of these 'arcs' are 
labelled in Fig. |4fb) with t heir morphological types as de- 
fined by lLaing etail (|2006bl) . In order to display them more 
clearly, Fig. [5] shows the relevant areas of the North and 
South jets on an expanded scale and with a narrower grey- 
scale range. The main characteristics of the arcs are as fol- 
lows. 



(i) The most prominent examples are at distances 
> 20 arcsec (7 kpc) from the nucleus, but there are hints of 
similar structures closer in. 

(ii) They all start from an edge of the jet and extend 
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Figure 5. Grey-scale images of 3C31 observed at 8440 MHz with a resolution of 0.75 arcsec FWHM. The grey-scale range is - 0.5 
beam -1 , (a) North jet, around the first bend; (b) South jet (note that the nucleus is not included in this panel). 
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Figure 7. Grey-scale of optical emission near 7000A observed 
with HST WFPC-2, superimposed on contours of total intensity 
at 8440 MHz with 0.25 arcsec FWHM resolution. The location of 
the dust dis k and the core radius of the inner hot gas c omponent 
are marked l|Martel et al.lll999l : Irlardcastle et alj|2002l ). 
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Figure 8. Vectors whose lengths are proportional to the mean de- 
gree of polarization, p, at L-band and whose position angles repre- 
sent the apparent magnetic field direction, x(0) + 7r/2, where x(0) 
is the zero- wavelength E- vector position angle from a rot ation- 
measure fit to all four L-band frequencies ( Laing et al J2008l) . The 
resolution is 5.5 arcsec FWHM and the polarization scale is indi- 
cated by the labelled bar. The diagram is designed to show the 
magnetic-field structure in the North tail over a region which is 
not sampled by the 4985- and 8440-MHz observations. 
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Table 1. Journal of observations. The columns are: (1) VLA 
configuration, (2) centre frequency, (3) bandwidth, (4) date of 
observation and (5) integration time. 



Conf v Au Date t 

MHz MHz min 



8460^ 


100.00 


1996 


Nov 12 


606 


1636 _ 


6.25 


1995 


Apr 28 


192 


1485 


6.25 


1995 


Apr 28 


192 


1435 


6.25 


1995 


Apr 28 


192 


1365 


6.25 


1995 


Apr 28 


192 


8440 


100.00 


1994 


Jun 6 


406 


8440 


100.00 


1994 


Jun 14 


412 


4985 


50.00 


1993 


Feb 25 


535 


1636 


25.00 


1994 


Jun 13 


192 


1485 


25.00 


1994 


Jun 13 


192 


1435 


25.00 


1994 


Jun 13 


192 


1365 


25.00 


1994 


Jun 13 


192 


8440 


100.00 


1994 


Dec 4 


242 


4985 


50.00 


1994 


Dec 4 


165 


1636 


25.00 


1994 


Nov 26 


140 


1485 


25.00 


1994 


Nov 26 


140 


1435 


25.00 


1994 


Nov 26 


139 


1365 


25.00 


1994 


Nov 26 


139 


8440 


100.00 


1995 


Apr 28 


69 


4985 


50.00 


1995 


Apr 28 


69 


1636 


25.00 


1995 


Apr 28 


61 


1485 


25.00 


1995 


Apr 28 


61 


1435 


25.00 


1995 


Apr 28 


69 


1365 


25.00 


1995 


Apr 28 


69 



a The slight difference in centre frequency used for the 
A-configuration observations has no measurable effect 
on the results and we refer throughout to the 
'8440-MHz observations'. 

" Most antennas were in the B configuration. 

inwards towards the jet axis in the downstream direction 
(i.e. away from the nucleus). 

(iii) The most prominent examples are narrow and ridge- 
like, but there are many fainter and broader features. 

(iv) Three of the arcs in the North jet clearly cross the 
axis. They are concave towards the nucleus and approxi- 
mately semi-circular i n form [the defining characteristics of 
type (i) arcs given bv lLaing et alj l|2006bl )]. 

(v) The remaining arcs in the North jet have higher 
brightness contrast at the outside edge of the bend. It is 
possible that they also extend across most of the jet, but at 
lower contrast. 

(vi) The features in the South jet have a greater degree of 
mirror symmetry across the jet and show similar brightness 
contrast on both sid es. They are type (ii) arcs as defined by 
lLaing et al.1 |2006b): oblique to the jet direction on either 
side of the jet but without large brightness gradients on- 
axis. They are located within the outer envelope of the jet 
emission, so the highest brightness gradients do not occur 
at the jet edges (Fig. [5p). 

(vii) The maximum local brightness enhancement is 



Table 2. Parameters of the images. 



V 


FWHM 


* 1 


^Jy beam 1 


MHz 


arcsec 


1 


QU 


8440 


0.25 


6 


6 


8440 


0.75 


7 


6 


8440 


1.50 


11 


6 


4985 


1.50 


12 


12 


1636 


1.50 


72 


64 


1485 


1.50 


73 


66 


1435 


1.50 


78 


68 


1365 


1.50 


76 


62 


Mean L-band 


1.50 


34 




8440 


5.50 


37 


12 


4985 


5.50 


26 


16 


1636 


5.50 


61 


22 


1485 


5.50 


52 


23 


1435 


5.50 


68 


25 


1365 


5.50 


68 


23 


Mean L-band 
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f»40% (for the brightest arc in the North jet); <20% for 
the other arcs. The implied emissivity enhancements in the 
arcs would be of order 400% to 200%, if their depths along 
the line of sight are of the same order as their apparent 
widths. 

We discuss the magnetic-field structure and other properties 
of the arcs in Section \4. 21 

Finally, Fig. [5] shows the inner jet region at 0.25- 
arcsec resolution in two different grey-scale representations. 
Fig. |6ja) shows that the brightest arc in the North jet is 
fully resolved, with a characteristic width «1 arcsec. The 
Figure also displays the knot at the base of the counter- 
jet which marks the brightening point on that side of the 
nucleus (we model the brightening as a sudden increase in 
emissivity at 2.5 arcsec from the nucleus in projection on 
both sides; lLaing fc Bridle! l2002ah . Fig. [IJb) illustrates the 
non-axisymmetric knot structure at the base of the North 
jet. 

3.4 Optical-radio superposition at high resolution 

HST WFPC - 2 obs ervations of 3C 31 were presented by 
iMartel et all (|l999l ). who showed that there is a complex 
central dust lane of radius ~3 arcsec. A superposition of 
our 8440-MHz image at 0.25 arcsec FWHM resolution on 
the HST image is shown in Fig. [7] In order to correct for un- 
certainties in the relative astrometry, we assumed that the 
brightest points on the optical and radio images coincide. 
The orientation an d size of the kpc-sc ale rotating molecu- 
lar disk studied by lOkuda et alj ((2005) are consistent with 
those of the dust lane. The radio jet flares and brightens 
on approximately the same (projected) scale as significant 
substructure (banding) in the dusty molecular disk (Fig. [7]), 
but the latter is grossly misaligned with the jet axis and 
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Figure 9. Vectors whose lengths are proportional to the mean 
degree of polarization, p, at 4985 MHz and whose position an- 
gles represent the apparent magnetic field direction, x(0) + 
7r/2, where x(0) is the zero-wavelength E-vector position an- 
gle from a rotation-meas ure fit to a ll frequencies between 1365 
and 4985 MHz (see lLaine et alj|2008t Fig.2c). The grey-scale is a 
Sobel-filtered, mean L-band image, normalized by total intensity. 
The resolution is 5.5arcsec FWHM and the polarization scale is 
indicated by the labelled bar. 



there is no evi dence for any direct interaction between the 
two structures l|de Kofi et alJ l200CO. The flaring scale is also 
comparable with t he core radius of t he inn er hot gas compo- 
nent identified by lHardcastle et alj (|2002h and indicated on 
Fig. [7] We have argued elsewhere that the flaring and recol- 
limation of the jets in 3C 31 is governed by the pressure of 
hot gas associated with the parent g alaxy, rather than by in - 
teraction with the cool component l|Laing fc Bridle! l2002bh . 



4 FIELD STRUCTURE AND INTENSITY 
GRADIENTS 

We have estimated the degree of polarization at zero fre- 
quency, p(0) and the direction of the apparent magnetic 
field (x(0)+7r/2, where x(0) is the zero- wavelength E-vector 
position angle) by correcting the observed values of p and 
the E-vector position angles for depolarization and Fara- 
day rotation, respectively. The depolarization in the North 
of the source at 5.5-arcsec resolution and over the entire 
area visible at higher resolution is very small and the change 
of E-vector position angle is accurately proportional to A 2 , 
so these corrections are well-determined (|Laing et al.| [2008). 
The rotation-measure distribution and its implications for 
the structur e of the mag n etoion ic medium around 3C 31 are 



analysed in lLaing et al.l (|200Sf ) ; here we describe only the 
apparent field structure within the jets and its relation to 
features of their total intensity. 



4.1 Field structure on large scales 

Fig. [8] shows the apparent magnetic-field structure of the 
outer North jet, spur and t ail derived from a n L-band RM 
fit at 5.5arcsec FWHM (see lLaing et al.ll2008l . Fig. 2b). The 
degree of polarization is high (p « 0.5 - 0.6) at the edges of 
the North tail and spur. There is a clear polarization min- 
imum in the centre of the tail, suggesting that the three- 
di mensional fie ld configuration might be close to model B 
of iLaingl ([1981), where the longitudinal and toroidal com- 
ponents are roughly equal but there is no field in the radial 
direction, perpendicular to the axis of the tail. 

Fig. \§\ shows the apparent magnetic-field structure at 
5.5-arcsec resolution from a 5-frequency fit. The remark- 
able straight edge of the South spur is highly polarized at 
4985 MHz, with an apparent field along the edge, suggesting 
that this part of the radio source has experienced compres- 
sion or shear at a large-scale planar boundary whose coun- 
terpart in the surrounding gas has yet to be identified. 

At 1.5-arcsec resolution (Fig. I10[) . the apparent field 
structure of the jets becomes clear. To a first approxima- 
tion, the outer regions of both jet s continue to show the 
configuration discussed in detail by lLaing fc Bridle] |2002a) 
for the inner 30 arcsec - an apparent field which is trans- 
verse on-axis but longitudinal at the edges of the jets. In the 
North jet (Fig. 110b ) this configuration persists for at least 
another 30 arcsec after the first bend. Thereafter, the on-axis 
field remains transverse, but there is insufficient sensitivity 
to measure the polarization at the jet edges. After the second 
bend, the apparent field becomes longitudinal over most of 
the jet. The degree of polarization is significantly enhanced 
at the outside edges of major bends in both jets. Where the 
South jet expands rapidly (Fig. HOfa ). it shows a predomi- 
nantly transverse apparent field, again with high p. There is 
a highly-polarized edge where the jet makes its final sharp 
turn prior to disrupting at the feature S in Fig. [2] 

The apparent magnetic fiel d in the brighte r jet at 0.25- 
arcsec resolution was shown by lLaing fc Bridle! l|2002al ). 



4.2 Magnetic fields in the arcs 

A superposition of the Sobel-filtered, 0.75-arcsec FWHM I 
image and vectors showing the magnitude of p and the ap- 
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Figure 10. Vectors whose len gths are proportional to p(0), the degree of polarization extrapolated to zero wavelength from a linear fit of p 
against A 2 IlLaing et al.l2008h . Their angles represent the a pparent magnetic fi eld direction, x(0) + 7r/2, where x(0) i s the zero-wavelength 
E-vector position angle corresponding to the RM image of lLaing et"afl 12008, Fig. 3c). The vectors are superposed on a grey-scale of the 
4985-MHz total intensity. The resolution is 1.5arcsec FWHM and the polarization scale is indicated by the labelled bars, (a) South jet; 
(b) North jet. 



parent field direction is shown in Fig. 1111 This reveals two 
striking effects: 

(i) The polarization is enhanced wherever there is a 
strong gradient in total intensity, particularly at the arcs. 
For example, p ~ 0.4 in the brightest arc in the North jet, 
even without correction for contamination by surrounding 
emission. 

(ii) The apparent field always appears parallel to the arcs. 

The arcs in 3C 296 sh ow similar, but l ess prominent 
magne tic-field struc t ures ( Laing et al.l l2006bh . 

In lLaing et al. I (|2006bh 7 we argued that the difference 
between type (i) arcs (found predominantly in the main jets) 
and type (ii) (usually in the counter-jets) was plausibly due 
to relativistic aberration. We suggested that a first-order 
model of an arc is a thin, axisymmetric shell of enhanced 
emissivity, concave towards the nucleus and travelling out- 
ward with approximately the velocity of the local mean flow. 
In that case, when our line of sight is tangential to part of 
the shell in its rest frame, we will see enhanced emission. 



Close to the centre-line of the jet, the shell can be approxi- 
mated as a planar sheet of material orthogonal to the axis in 
the rest frame of the flow. Relativistic aberration causes the 
sheets to appear differently orientated to the line of sight in 
the main and counter-jets so that those in the main jet can 
be observed nearly edge-on while those in the counter-jet 
are always significantly rotated about a line perpendicular 
to the jet axis in the plane of the sky and are therefore less 
prominent both in intensity and in brightness gradient where 
they cross the jet axis. These differences between the main 
jet and counter-jet are less pronounced for the slower flow 
at the edges of the jets, where the arcs can therefore appear 
more symmetrical - as observed. If the magnetic field in the 
arcs is compressed to lie in the plane of the shells, we also 
expect to see a high degree of polarization with the app arent 
field direction along the arc, as observed i|Laing||l980l ). 



The arcs therefore provide further support for the hy- 
pothesis that the main differences in appearance between 
the inner main jet and counter-jet in 3C 31 arise from rel- 



12 R.A. Laing et al. 




01 07 27.5 27.0 2 6.5 26J0 25.5 25.0 01 07 26.0 25,5 25.0 24.5 24.0 25.5 

RIGHT ASCENSION (J2000) RIGHT ASCENSION (J2Q00) 

Figure 11. The colour plots show Sobel-filtered, 8440-MHz / images at a resolution of 0.75arcsec FWHM superimposed on vectors 
whose magnitudes are proportional to the degree of polarizatio n, p, and whose di rections are those of the apparent magnetic field. The 
field directions were derived by interpolating the RM image of lLaing et al. | teOOSl. Fig. 3c) onto a finer grid and using it to correct the 
observed 8440-MHz E-vector position angles to zero wavelength, (a) South jet; (b) North jet. 



ativistic bulk motion in intrinsically similar jet outflows on 
the two sides of the nucleus. 



5 SPECTRA 

5.1 Comparison of images at different frequencies 

At 4985 and 8440 MHz, we cannot observe the full range 
of spatial frequencies covered by the D configuration at L 
band. We have therefore investigated the sensitivity of the 
spectral index determination between 1365 and 4985 MHz 
at 5.5 arcsec resolution to differences in the (u, v) coverage 
by comparing the spectral index distributions obtained from 
images made with: 

(i) all of the observed data at both frequencies, imaged 
with zero-spacing flux densities specified; 

(ii) matched coverage at the two frequencies, with the 
central portion of the (u, v) plane, including the zero spac- 
ing, deleted at 1365 MHz); 

(iii) as in (ii), but with the zero-spacing flux densities 
included. 



We show results only for areas of 3C 31 where the spectral 
indices obtained by these three methods agree (Ao < 0.1). 
The limitations in (u, v) coverage and primary beam size im- 
ply that we can derive reliable spectra only for the North and 
South jets and the South spur between 1365 and 4985 MHz 
at a resolution of 5.5 arcsec: the tails and the North spur 
are seen only at L-band, within which the range of observed 
frequencies is too small for us to derive reliable spectral in- 
dices. 

Similar considerations apply at 1.5-arcsec resolution, 
where we have also compared spectral-index distributions 
derived from CLEAN and maximum-entropy deconvolutions. 
Spectra derived from the 8440-MHz images are reliable only 
within ±70 arcsec of the nucleus, for the following reasons: 

(i) As mentioned earlier, the short-spacing coverage of the 
VLA D-config uration is inadequate to image structure on 
larger scales ^Taylor. Ulvestad fc Perlevll2003r ). 

(ii) A comparison of spectral indices derived from im- 
ages at 5.5 and 1.5-arcsec resolution, and from CLEAN and 
maximum-entropy images at the higher resolution, showed 
differences at the level of Aq > 0.1 farther from the nucleus. 

(iii) In order to keep the rms error in fitted spectral in- 
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Figure 12. (a) Sobel-filtered, mean L-band image (Section 12.21 . normalized by total intensity, at a resolution of 5.5arcsec. (b) and (c) 
Spectral index, a from weighted least-squares, power-law fits to the total intensity, (b) 5-frequency fit between 1365 and 4985 MHz at a 
resolution of 5.5 arcsec FWHM. Only points with / > 5<x at all frequencies are shown. The colour range is 0.5 < a < 1.66. (c) 6-frequency 
fit to the maximum-entropy images at a resolution of 1.5 arcsec FWHM for the inset area in panel (b). Only points with / > 10(7 at all 
frequencies were used in the fit. The colour range is 0.5 < a < 0.8. The boxes mark the regions used for the fits in Fig. 1 13 1 and Table|3] 
All of the images have been rotated anticlockwise by 19.7 deg. 



dex <0.1 at 1.5-arcsec resolution, we required that all input 
images had / > 10a; this condition was satisfied only in the 
central regions. 

For the same reasons, we can search for deviations from a 
power-law spectrum only in the inner jets. 

Experience suggests that the combination of calibration 
and deconvolution errors results in an rms oi ~ 0.03/. When 
fitting spectra, we have therefore added this term and the 
off-source noise from Table [2] (corrected for primary-beam 
attenuation) in quadrature. As noted in Section 12.21 chro- 
matic aberration does not affect the measured brightnesses 
for 3C 31 significantly. The reason is that the effect preserves 
flux density and 3C 31 is very well resolved on all scales 
where it is important. 



5.2 Power-law spectra 

The best representation of the spectrum at a resolution 
of 5.5 arcsec is given by a weighted least-squares fit of 
power-law spectra to the / images over the frequency 
range 1365 - 4985 M Hz (Fig. fT^bt. The 'spectra l to- 
mography' technique o f iKatz-Stone fc Rudnickl (|l997T ) and 
iKatz-Stone et ahl (1 19991 ) leads to similar conclusions, but our 
use of all of the available frequencies minimizes noise. 

The spectral index variation over the inner jets is shown 
in more detail in Fig. I12f cl . which is derived from a six- 



frequency power-law fit at 1.5-arcsec resolution. Points are 
only plotted if: (a) / > 10a at all frequencies (resulting in 
a maximum rms error in spectral index of 0.11) and (b) the 
spectral indices derived from CLEAN and maximum-entropy 
images differ by <0.075. In order to search for deviations 
from power-law spectra at high signal-to-noise ratio, we have 
integrated the flux densities over the rectangular regions 
shown in Fig. I12f c) for the 1.5-arcsec FWHM images, using 
the same blanking criterion as for the spectral-index map in 
Fig. 112( 0). The results are shown in Fig |13l and the spectral 
indices are tabulated in Table [3] 



5.3 Spectral variations 

5. 3. 1 Description 

The principal features of the spectra shown by Figs [12] and 
[T3]are as follows: 

(i) There is no evidence for any significant deviation from 
power-law spectra wherever we have reliable 8440-MHz data 
(i.e. within 70 arcsec of the nucleus). All of the fits shown 
in Fig. 1131 except for SI have a probability > 0.45 that the 
chi-squared statistic would exceed the measured value by 
chance, suggesting that we have slightly overestimated the 
multiplicative error term. For SI (the faint region of the 
counter-jet close to the core) there are significant residual 
artefacts at L band which are not represented in our noise 
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Figure 13. Power-law fits to the total-intensity images over the 
rectangular regions shown in Fig. H2f cV (a) North jet, (b) South 
jet. 



estimate. Our constraints on spectral curvature are limited 
by calibration uncertainties at the 1 - 2% level (it is clear 
from Fig.[[3]that the 1636- and 4985-MHz points are slightly 
low relative to the best-fitting power law and the 8440-MHz 
points slightly high). In any case, the curvature expected 
over our restricted frequency range is small, because of the 
finite width of the synchrotron emission spectrum for a sin- 
gle electron. For these reasons, we will not discuss spectral 
curvature in what follows. As the calibration differences are 
common to all of the fits, the differences in spectral index 
between regions are better determined than the absolute val- 
ues or curvature, however, with an estimated rms error of 
only 0.01. 

(ii) Within «7 arcsec of the nucleus in both jets, the 
spectral index at 1.5-arcsec resolution is slightly steeper 
((a) = 0.62) than the average for the inner jets. 

(iii) From 7 - 50arcsec in both jets, the mean spectral 



index is in the range 0.55 - 0.57. Farther from the nucleus, 
there is a gradual spectral steepening. 

(iv) There are also slight, but significant variations in 
spectral index across both jets within ~30 arcsec of the nu- 
cleus in the sense that their West edges tend to have flatter 
spectra ((a) = 0.52 - 0.54; Fig. [12}:). We are confident that 
this is not an instrumental effect: it is clearly visible at both 
resolutions and on the CLEAN and maximum-entropy images. 
Any residual misalignments of the images, particularly at 
1.5-arcsec resolution, should be extremely small. The most 
likely remaining error (confusion with the unresolved radio 
core at the base of the North jet) would not produce a dif- 
ference in spectral index across the jets. Flatter spectra are 
also se en on both edges of the jets in NGC315 l|Laing et al.l 
l2006al ). 

(v) The spectra of the two brightest arcs in the North jet 
are also slightly flatter than that of the surrounding emission 
(Fig.[T2j;; Table |3j). 

(vi) There is a clear spectral separation between the 
South jet and the spur emission, matching the separa- 
tion of these regions defined by the sharpest brightness 
gradients (Figs 112b and b). This is particularly clear 
where the jet first enters the diffuse emission. The spec- 
tral identity of the jet is evidently maintained even af- 
ter it bends abruptly about 2 arcmin South of the nu- 
cleus and remains until it terminates in a region of high 
brightness gradient. The outer South jet in 3C 31 is there- 
fore a clear example of th e type of spectral structure 
noted in other FRI sources (|Katz-Stone fe Rudnicklll997i ; 
iKatz-Stone et all 1 19991 ) wherein a flatter-spectrum 'jet' 
with a distinct spectral identity is superposed on steeper- 
spectrum 'sheath' emission. 

(vii) There are hints of a similar spectral difference be- 
tween the jet and the spur and tail in the North, but the 
sensitivity and short-spacing coverage of the 4985-MHz ob- 
servations are inadequate to confirm this. The difference is 
seen clearly in WSRT observations at 0.6 and 1 .4 GHz with a 
resolution of 29 x 55 arcsec FWHM (fig. 3g of |jagerdfl987l) . 

(viii) Approximate subtraction of the steeper-spectrum 
sheath emission by linear interpolation under the jets shows 
that there must be some intrinsic steepening of the jet spec- 
trum with increasing distance from the nucleus in addition 
to the apparent steepening caused by emission from sheath 
material along the line of sight. The complex morphology 
of the South spur and tail make it difficult to separate the 
spectral components unambiguously, however. 

(ix) At 5.5-arcsec resolution, the spectral index of the 
North jet is flatter by 0.1 - 0.2 on-axis than at the edges 
in the region of the double bend, between 80 and 120 arcsec 
from the nucleus, the flattest spectrum being associated with 
a bright filament (Fig. 112b ). Similarly, transverse spectral 
gradients in the opposite sense to those described in (iv) - 
i.e. with the jet spectrum flatter on-axis - are seen at dis- 
tances > 50 arcsec from the nucleus at 1.5-arcsec resolution 
(Fig. 112b ). These effects are most likely to be produced by 
blending of flat-spectrum emission from the jets with steep- 
spectrum emission from the surrounding spurs and tails. 

(x) lAndernach et alj <| 19921 ) found that within 9 arcmin of 
the core the spectrum observed with lo w resolution fl attens 
above 5 GHz, as previously suspected bv lBurcb] (1 19771 ). They 
interpreted this effect as a superposition of components with 
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Table 3. Spectral indices from 6-frequency power-law fits for the 
regions shown in Fig, ll2f c), The flux densities used in the fits were 
derived by integration over the images at 1.5-arcsec resolution. 
The rms errors on the spectral indices are all between 0.017 and 
0.020, and are dominated by calibration uncertainties. The rms 
errors on spectral-index differences between regions are faO.01. 



Region 

N arc 2 

N arc 1 

N edge 

N4 

N3 

N2 

Nl 

SI 

S2 

S3 

S4 

S5 

S edge 



a 

0.54 
0.56 
0.52 
0.56 
0.55 
0.58 
0.62 
0.62 
0.58 
0.57 
0.57 
0.63 
0.54 



different spectra, a conclusion which is confirmed directly by 
our data. 



5.3.2 Spectral variations in FRI jets 

These observations contribute to the developing picture of 
spectral variations in the bases of FRI jets, as follows. 

(i) Where the jets first brighten, there is a remarkably 
small dispersion around a spectral index of a — 0.62 in 
the three sources w e have studie d in detail: 3C 31, N GC 315 
i|Laing et alj|2006al) and 3C 296 (lLaing et alj[2006bh . as well 
as 3C 66B (|Hardcastle et ail l200lf ) . The average is domi- 
nated by emission immediately after the brightening point 
(Fig.®. 

(ii) The spectral index of the fainter emission close 
to the nucleus in 3C 449 dRudnick fc Katz-Stond 
I1997D, PKS1333-33 jKilleen et all Il986h and 3C66B 
1 Hardcastle et al.l l200lf ) appears to be comparable or 
slightly steeper although the uncertainties are larger. We 
do not have adequate resolution to measure the spectral 
index in these parts of the jets in 3C 31. 

(iii) Farther from the nucleus, the radio spectra flatten 
slightly to a = 0.50 - 0.55. This flattenin g occurs where 
our kinematic models require deceleration |Laing fc Bridle! 
l2002al ; ICanvin et al1l2005l ; lLaing et al.ll2006bl) . 

(iv) A related result is that an asymptotic low-frequency 
spectral index of 0.55 is co mmon in FRI jets over larger 
areas than we consider here ll Young et al.ll2005h . 

(v) Flatter-spectrum edges can be seen where the jets 
are isolated from signi ficant surrounding diffuse emission. 
Our kinematic models |Laing fc Bridlell2002al ; ICanvin et al.l 
120051 ; lLaing et al. 2006b) show that all of the jets have sub- 
stantial transverse velocity gradients in these regions. 

(vi) Our kinematic models and spectral measurements to- 
gether suggest a relation between spectral index and flow 
speed f3 = v/c, ranging from a = 0.62 for fa 0.8 to a m 0.5 
for < 0.2. Spectra would then be expected to flatten with 



distance from the nucleus (as the jets decelerate) and from 
centre to edge (as a result of transverse velocity gradients). 

(vii) Synchrotron X-ray emission from the main jets in 
3C66B, 3C31, 3C 296 and NGC 315 is strongest relative to 
the radio close to the nucleus, at or before the brighten- 
ing p oint l|Hardcastle et al.ll200ll , I2OO21 . 120051 ; IWorrall et ail 
2007). X-ray emission is still detected from the flatter- 
spectrum regions farther out, but at a lower level relative 
to the radio. 

(viii) There is approximate morphological correspondence 
between the radio and X-ray brightness distributions, but 
they differ in detai l. In the best-resolved case, NGC 315 
jWorrall et alj |2007). the X-ray emission is clearly extended 
across the jet. Particle acceleration appears to be distributed 
throughout the jet volume, rather than being exclusively as- 
sociated with discrete knots or with the boundary. 

(ix) The ratio of X-ray to radio emission decreases and 
the radio spectrum starts to flatten where our kinematic 
models show that the jets decelerate rapidly from speeds of 
w 0.8 - 0.9 jLaing fc Bridldl2002al; lHardcastle et al.ll2o"oi 
ICanvin et aT]|2005l ; IWorrall et al.ll2007l ). 

(x) The strongest X-ray emission is not associated with 
the flattest radio spectra, but rather with some particle 
acceleration process whose characteristic energy index is 
s — 2a + 1 = 2.24. This is intriguingly close to the asymp- 
totic value of s = 2.23 for first-order Fermi acceleration 
at relativistic shocks in the limit of large shock (or, equiv- 
alently, upstream) Lorentz facto r V, a result whi ch holds 
even for relativi stical l y hot jets (|Kirk et al.l |2000| and ref- 
erences therein; iKirkl |2005h . The asympt otic value of s is 
only approached for T = ( 1 -/? 2 )~ 1/2 > 10 (|Kirk et al.ll2000l ; 
iLemoine fc Pelletierl [2003) . however, so we would have to 
suppose that some fraction of the flow before the brighten- 
ing point has a bulk Lorentz factor F > 10 and decelerates 
to r ~ 2, as required by our kinematic models. The shock 
would then have to be oblique, since the normal compo- 
nent of the the incoming plasma speed would be reduced 
to ps 1/3 at the shock. Such a flow would be difficult to 
detect because of the narrow beaming angle, even if it con- 
tains relativistic particles. Emission before the brightening 
point would then come p rimarily from a s low su rface layer, 
as suggested for 3C 31 by lLaing fc Bridld (|2002al ). 

(xi) As the jet slows down and the shocks become both 
less relativistic and weaker, the energy index would decrease, 
eventually reaching s = 2 (a = 0.5) when the shocks be- 
come non-relativistic (|BellHl978h , Although this provides a 
natural reason for the spectrum of a jet to flatten as it de- 
celerates, the required flow speeds are higher than we infer. 
In particular, the spectral index of a — 0.61 is seen over a 
distance of ~8 kpc (corrected for projection) along the jet 
axis in NGC 315, i n a region where we i nfer the on-axis flow 
to have V « 2.1 jCanvin et all 120051 ; lLaing et al] l2006a[ ). 
We would then have to suppose that some part of the flow 
remains ultrarelativistic (and effectively invisible) on large 
scales. 

(xii) A second possibility is that the process t hat produces 
the fla tter spectrum is associated with shear l|Laing et al.l 
l2006al ). In 3C31, flatter-spectrum edges occur predomi- 
nantly on the Western edges of the jets, i.e. on the outer 
edges of bends, perhaps consistent with this idea. Second- 
order Fermi acceleration driven b y turbulence in the shear 
layer (|Stawarz fc Ostrowskil I2002T ) appears to be a viable 
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process in FR I jets, but the s hear a cceleration mechanism 
described by iRieeer fc Duffvl (|2004) is unlikely to be effi- 
cient enough to accelerate the X-ray-emitting electrons in 
the relatively modest velocity gradients we infer. 

(xiii) The idea that two different acceleration processes 
are required has also been suggested on the basis of evidence 
from the X-r ay morphology and spectrum of the brighte r 
jet in CenA (jKataoka et all l200rj ; lHardcastle et al.l l2007t) . 
Here, the compact knots observed close to the nucleus are 
thought to be associated with shocks, while a truly diffuse 
acceleration mechanism dominates at larger distances. 

Far from the nucleus, the picture is complicated by the 
presence of diffuse emission surrounding the jets. As well 
as a smooth steepening of the jet spectrum (expected from 
synchrotron and adiabatic losses affecting a homogeneous 
electron population), multiple spectral components are ob- 
served. Jets appear to retain their identities even after en- 
tering regions of diffuse emission and are clearly identifi- 
able by their flatter spectra. They are usually separated 
from the surrounding emission by sharp brightness gradi- 
ents. This separation is observed in FR 1 sources with bridges 
of emission extendi ng back towards the nucleus (e.g. 3C 296; 
lLaing et al.ll2006bT ) as well as sources with tails and spurs 
like 3C3f. Although there is an overall trend for the spec- 
trum of the diffuse emission to stee pen towards the nu cleus 
in bridges and away from it in tails (jParma et al.l[f999l ). the 
variations in individual objects such as 3C 3f are complex. 
The termination regions of jets in tailed FRI sources are 
perhaps best regarded as bubbles which are continually fed 
with fresh relativistic plasma by the jets and which in turn 
leak material into the tails. Their spectral steepening would 
then be governed by a combination of continuous injection, 
adiabatic, synchrotron and inverse Compton energy losses 
and escape. 



6 SUMMARY AND CONCLUSIONS 

Our new images of the FRI radio galaxy 3C3I reveal a 
complex substructure of arcs and filaments within the jets 
and tails that is clearly related to the organization of the 
apparent magnetic field structure. There is a strong corre- 
spondence between steep brightness gradients and the or- 
ganization of the apparent field direction perpendicular to 
suc h gradients. 

lLaing fc Bridie] (|2002al ) have shown that the detailed 
brightness and polarization properties of the inner jets in 
3C 31 are consistent with interpreting them as intrinsically 
symmetrical relativistic flows that decelerate with increasing 
distance from the nucleus, so that the observed asymmetries 
also decline with distance. We suggest that the observed 
differences between the intensity and polarization structures 
of the arcs in the inner main and counter-jets are also caused 
by the effects of aberration on emission from intrinsically 
similar structures. 

Our new data further suggest that we are observing the 
transition between predominantly relativistic asymmetries 
which dominate in the straight inner parts of the main jet 
and the counter-jet and environmental asymmetries which 
dictate the appearance of 3C 31 at distances > 1 arcmin from 
the nucleus. Effects of an asymmetric environment include: 
(a) the different scales on which bending begins in the North 



and South jets; (b) the dissimilar morphologies of the North 
and South spurs and tails, and (c) the lack of any systematic 
correspondence between the locations, shapes and bright- 
nesses of substructures in the outer parts of the North and 
South jets. 

Our images also show that both jets in 3C 31 retain their 
spectral identity far from the galactic nucleus, where their 
emission is superimposed on larger-scale, steeper- spectrum 
structures. Spectral imaging and intensity-gradient imaging 
(Sobel filtering) techniques both delineate the same bound- 
aries between the jets and their environs on both the North 
and South sides of the galaxy. 

The radio spectra of the inner jets of 3C 31 and other 
FRI radio galaxies exhibit systematic trends which are re- 
lated both to jet kinematics and to X-ray emission: 

(i) The characteristic radio spectral index where the jets 
first brighten is a = 0.62 ± 0.01 . The ratio of X-ray/radio 
emission is highest at and before the brightening point. 

(ii) Farther out, the spectrum flattens slightly, to a ~ 
0.50 — 0.55 and the X-ray emission is fainter relative to the 
radio. The onset of this spectral flattening coincides with 
deceleration of the flow in our kinematic models. 

(iii) Flatter spectra are seen at the jet edges where the 
jets are isolated from significant diffuse emission, in regions 
where our kinematic models require significant transverse 
velocity gradients. 

The correlation between radio spectral index and X-ray 
emission suggests a direct association with the particle- 
acceleration process. The energy index of 2.24 corresponding 
to a = 0.62 is close to the asymptotic value for relativistic 
shocks in the limit of high Lorentz factors, and the spec- 
tral flattening could result from progressive weakening of 
the shocks as the jets slow down. The Lorentz factors of at 
least some parts of the flow would have to be much higher 
than we infer, however. Alternatively, there may be two dif- 
ferent acceleration mechanisms operating on kpc scales in 
these jets, one dominant from the nucleus until the region 
where jets brighten abruptly in the radio, the other taking 
over at larger distances and perhaps associated with trans- 
verse velocity shear. 



ACKNOWLEDGEMENTS 

RAL would like to thank the Istituto di Radioastronomia, 
NRAO and Alan and Mary Bridle for hospitality during the 
course of this work. We acknowledge travel support from 
NATO Grant CRG93I498. We thank John Kirk for advice 
on particle acceleration mechanisms and the referee, John 
Wardle, for his helpful report. The National Radio Astron- 
omy Observatory is a facility of the National Science Foun- 
dation operated under cooperative agreement by Associated 
Universities, Inc. 



REFERENCES 

Andernach, H., Feretti, L., Giovannini, G, Klein, U., Ros- 

setti, E., Schnaubelt, J., 1992, A&AS, 93, 331 
Arp, H., 1966, ApJS, 14, I 
Bell, A.R., 1978, MNRAS, 182, 147 



Observations of 3C 31 17 



Burch, S.F., 1977, MNRAS, 181, 599 
Burch, S.F., 1979, MNRAS, 187, 187 

Butcher, H.R., van Breugel, W.J.M., Miley, G.K., 1980, 
ApJ, 235, 749 

Canvin, J.R., Laing, R.A., Bridle, A.H., Cotton, W.D., 

2005, MNRAS, 363, 1223 
Croston, J.H., Birkinshaw, M., Conway, E., Davies, R.L., 

2003, MNRAS, 339, 82 
de Koff, S., Best, P., Baum, S. A., Sparks, W., Rottgering, 

H., Miley, C, Golombek, D., Macchetto, F., Martel, A., 

2000, ApJS, 129, 33 
De Vaucouleurs, C, De Vaucouleurs, A., Corwin, H.G. Jr., 

Buta, R., Paturel, G., Fouque, P., 1991, Third Reference 

Catalogue of Bright Galaxies, Springer- Verlag, New York 
Fanaroff, B.L., Riley, J.M., 1974, MNRAS, 167, 31P 
Fomalont, E. B., Bridle, A. H., Willis, A.G., Perley, R. A., 

1980, ApJ, 237, 418 
Fraix-Burnet, D., Golombek, D., Macchetto, F., 1991, AJ, 

102, 562 

Hardcastle, M.J., Birkinshaw, M., Worrall, D.M., 2001, 

MNRAS, 326 1499 
Hardcastle, M.J., Worrall, D.M., Birkinshaw, M., Laing, 

R.A., Bridle, A.H., 2002, MNRAS, 334, 182 
Hardcastle M.J., Worrall D.M., Birkinshaw M., Laing R.A., 

Bridle A.H., 2005, MNRAS, 358, 843 
Hardcastle M.J., et al., 2007, ApJ, 670, L81 
Huchra, J.P, Vogeley, M.S., Geller, M.J., 1999, ApJS, 121, 

287 

Jagers, W.J.M., 1987, A&AS, 71, 75 

Kataoka, J., Stawarz, L., Aharonian, F., Takahara, F., Os- 
trowski, M., Edwards, P.G., 2006, ApJ, 641, 158 

Katz-Stone, D.M., Rudnick, L., 1997, ApJ, 488, 146 

Katz-Stone, D.M., Rudnick, L., Butenhoff, C, 
O'Donoghue, A. A., 1999, ApJ, 516, 716 

Keel, W.C., 1988, ApJ, 329, 532 

Killeen N.E.B., Bicknell G.V., Ekers R.D., 1986, ApJ, 302, 
306 

Kirk, J.G., 2005, in Chen P., Bloom E., Made- 
jski G., eds, 22nd Texas Symposium on Rela- 
tivistic Astrophysics, at Stanford. SLAC, Stanford, 
|http: / /www. slac.stanford.edu/econf/C041213 / 1 

Kirk, J.G., Guthman, A.W., Gallant, Y.A., Achterberg, A., 
2000, ApJ, 542, 235 

Komossa, S., Bohringer, H., 1999, A&A, 344, 755 

Laing, R.A., 1980, MNRAS, 193, 439 

Laing, R.A., 1981, ApJ, 248, 87 

Laing, R.A., Bridle, A.H., 2002a, MNRAS, 336, 328 

Laing, R.A., Bridle, A.H., 2002b, MNRAS, 336, 1161 

Laing, R.A., Bridle, A.H., 2004, MNRAS, 348, 1459 

Laing, R.A., Bridle, A.H., Parma, P., Murgia, M., 2008, 
MNRAS, submitted 

Laing, R.A., Canvin, J.R, Cotton, W.D., Bridle, A.H., 
2006a, MNRAS, 368, 48 

Laing, R.A., Canvin, J.R, Bridle, A.H. Hardcastle, M.J., 
2006b, MNRAS, 372, 510 

Lara, L., Cotton, W.D., Feretti, L., Giovannini, G., Venturi, 
T., Marcaide, J.M., 1997, ApJ, 474, 179 

Lemoine, M., Pelletier, G., 2003, ApJ, 589, L73 

Martel, A., Baum, S. A., Sparks, W. B., Wyckoff, E., 
Biretta, J. A., Golombek, D., Macchetto, F. D., de Koff, 
S., McCarthy, P. J., Miley, G. K., 1999, ApJS, 122, 81 

Okuda, T., Kohno, K., Iguchi, S., Nakanishi, K., 2005, ApJ, 



620, 673 

Parma, P., Murgia, M., Morganti, R., Capetti, A., de 

Ruiter, H.R., Fanti, R., 1999, A&A, 344, 7 
Pratt, W.K., 1991, Digital Image Processing, Wiley, New 

York 

Rieger F.M., Duffy P., 2004, ApJ, 617, 155 
Rudnick, L., Katz-Stone, D.M., 1997, ApJ, 488, 146 
Smith, R.J., Lucey, J.R., Hudson, M.J., Schlegel, D.J., 

Davies, R.L., 2000, MNRAS, 313, 469 
Stawarz, L., Ostrowski, M., 2002, ApJ, 578, 763 
Strom, R.G., Fanti, R., Parma, P., Ekers, R.D., 1983, A&A, 

122, 305 

Taylor, G. B., Ulvestad, J.S., Perley, R.A., The Very Large 
Array Observational Status Summary, NRAO, Socorro, 
URL: 

http: //www. vla.nrao . edu/astro/guides/vlas/ current/ 
Wardle, J.F.C., Kronberg, P.P., 1974, ApJ, 194, 249 
Worrall D.M., Birkinshaw M., Laing, R.A., Cotton, W.D., 

Bridle, A.H., 2007, MNRAS, 380, 2 
Young A., Rudnick L., Katz D., Delaney T., Kassim N.E., 
Makishima K., 2005, ApJ, 626, 748 
Zwicky, F., Kowal. C. T., 1968, Catalogue of Galaxies and 
Clusters of Galaxies, Vol.6, California Institute of Tech- 
nology, Pasadena 



